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ABSTRACT 

We consider axisymmetric relativistic jets with a toroidal magnetic field and an 
ultrarelativistic equation of state, with the goal of studying the lateral structure of jets 
whose pressure is matched to the pressure of the medium through which they propagate. 
We find all self-similar steady-state solutions of the relativistic MHD equations for this 
setup. One of the solutions is the case of a parabolic jet being accelerated by the pressure 
gradient as it propagates through a medium with pressure declining as p{z) oc z~'^. As 
the jet material expands due to internal pressure gradients, it runs into the ambient 
medium resulting in a pile-up of material along the jet boundary. The magnetic field 
acts to produce a magnetic pinch along the axis of the jet. The relative importance 
of the two components - the boundary pile-up and the magnetic pinch - depends on 
the magnetization parameter of the jet and on the total energy it carries relative to 
the ambient pressure. Such jets can be in a lateral pressure equilibrium only if their 
opening angle 9j at distance z is smaller than about I/7, where 7 is the characteristic 
bulk Lorentz-factor at this distance; otherwise, different parts of the jet cannot maintain 
causal contact. We construct maps of optically thin synchrotron emission from our self- 
similar models to illustrate the observational features of such jets. We suggest that the 
boundary pile-up may be the reason for the limb-brightening of the sub-parsec jet of 
M87. We find that if the synchrotron emissivity falls with the distance from the jet axis 
(as in the magnetically pinched case), the polarization fraction rises toward the edge, 
as seen in 3C273 and Mkn501. Projection effects and the emissivity pattern of the jet 
have a strong effect on the observed polarization signal, so the interpretation of the 
polarization data in terms of the geometry of magnetic fields is rather uncertain. For 
example, jets with toroidal magnetic fields display the 'spine-sheath' polarization angle 
pattern seen in some BL Lac objects. 
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Introduction 



Accretion onto black holes is frequently accompanied by collimated outflows. Although many 
of the physical properties of these flows proved difficult to measure, strong evidence exists that 
they are relativistic in a variety of environments. This evidence includes the one-sided appear- 
ance of jets and observatio ns of apparent superluminal motions in Galactic black-hole binaries 



(jMirabel &: Rodrigueall999l ) and in active g alactic nuclei (AGN: e.g., 
as the rapid gamma-ray variability in AGN (|Dondi &: Ghisellinilll995l : 



Jorstad et al 



Aharonian et al 



2Qm). a s well 



20051 ^. 



Jets are also magnetized. The polarization fraction of s ynchrotron radiation from AGN jets ap- 
proaches its theoretical maximum value in some objects (e.g., Zavala Sz TavloijboOS : Lister Sz Homan 



2005 



Pushkarev et al.ll2005l), indicating that t he magnetic field is h ighly ordered. The maps of 



rotation measure of 3G273 (Asada et al. 



2002 



Zavala &: Tavloii l2005l ) and of several BL Lac ob- 



jects (jGabuzda et al.l 12004 ) show strong gradients of the rotation measure across the jet which 
are expected if the toroidal component dominates. The polarization position angles indicate 
toroidal fields in more than half of all BL Lac objects (jGabuzda et al.ll2000l : iLister &: Homanll2005l : 
Pushkarev et al.ll2005l). but such observations may be subject to projection effects and relativistic 
aberration (jPariev et al.l 120031 : iLyutikov et al.ll2005l ). The strength of the magnetic field is poorly 
constrained, so its true dynamical importance cannot be evaluated from observations. 

The current theoretical paradi gm is that most relativist i c jets are pro b ably accelerated b y 
large-scale electromagnetic stresses (jBlandford &: Paynd Il982l : iLi et al.l Il992l : IContopoulod Il994l ) . 
Acceleration is rapid until the flow passes the fast magnetosonic surface and attains rough equipar- 
tition between kinetic energy flux and Poynting flux. Further increase in the jet Lorentz-factor 
tends to be more gradual, but eventually it is p ossible for most of the Po yn ting flux to be con- 
verted to kinetic energy (IBegelman &: Lil Il994l : IVlahakis &: Konigll l2003al lbl: iBeskin &: Nokhrina 



2006: Komissarov et al 



20071 ). Ultimately, the flow is confined laterally by either the inertia of 



the disk at the base of the jet, in which case the collimating force is tra nsmitted along t he jet by 
magnetic tension, or the pressure of the ambient medium along its walls (jBegelmanl Il995l ) . 



In this paper we investigate the interior of laminar ultrarelativistic jets whose boundary pres- 
sures are matched to their surroundings. The flux-freezing condition implies that the poloidal 
component of the mag netic field decays faster than the toroidal component in an expanding jet 
([Begelman et al.lll984l ). and we consider the structu re of the jet f a r eno ugh from its initial accel- 
eration site that the toroidal component dominates. IContopoulosI (|l995l ) showed formally that at 
large distances down the jet the poloidal magnetic field is dynamically unimportant and discussed 
the critical points of such solutions. 

As forces exterior to the jet gradually collimate it, curving the outermost streamlines slightly 
toward the axis, some of the material in the jet interior runs into the wall and piles up. Thus, jets 
may develop an edge-brightened or 'hollow-cone' structure, with much of the jet energy, momen- 
tum and emissivity concentrated into a narrow sheath along the interface between the jet and its 
surroundings. The pinching effect due to magnetic tension has the opposite effect, tending to focus 
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jet en ergy toward the axis (jUstvugova et al.l Il995l . 120001 : IStone k. Harded I2OOOI : iKomissarov et al 



20071 ). We find that for some functional dependencies of the ambient medium pressure Pa{z) such 
jets can develop self-similar structures that admit a wide range of prescribed entropy and toroidal 
magnetic field distributions, and we investigate the interplay of pressure and magnetic field for the 
self-similar solutions. 

In Section [2] we describe the basic equations, our assumptions and the family of self-similar 
solutions that we find. In Section [3] we discuss the appearance of the synchrotron radiation from 
our solutions. We conclude in Section HI In Section [21 unless otherwise noted we use units in which 
c = 1 and the magnetic field is redefined to eliminate factors of 47r from Maxwell's equations. The 
magnetic field in Gaussian units is times our magnetic field B and the current (and its density) 
in Gaussian units is 1/ \/47r times our current (and density) . Our Greek indices run from to 3 and 
the Latin ones from 1 to 3, and summation is assumed over repeating indices (using the Minkowski 
metric for Greek indices). We use Gaussian and CGS units in Section [31 



Self-similar solutions 



2.1. Main equations 



We start from the covariant form of the equations of r elativistic MHD in flat space-time (jPixon 



1978: 


Anile 


1989: 


Komissarov 


1999: 



which the pressure tensor is isotropic: 



Vuipu"") = 0; 
V^T^^ = 0; 
V^(6^'u^ - V'u^') = 0. 



(1) 
(2) 
(3) 



These equations are, from top to bottom, the continuity equation, the energy and momentum 
conservation equations, and the source-free part of Maxwell's equations. Here, p is the mass 
density, = (7, V7) is the 4-velocity of a fluid element. 



T^'^ = {^p + £ + P + b^)u^u'' + {P + h^)g^'' 



(4) 



6^6^ is the 



is the energy tensor, e is the internal energy density per unit volume, P is pressure, 6^ 
electromagnetic energy density and g^'^ is the space metric (flat in our case). If B is the magnetic 
field (in the observer's frame) of the fluid element moving with u^, the components of the 4- vector 
bf^ are W = B''"u'' and U = (i?* + b^u^)/u^, so the source-free Maxwell's equations ([3]) are the same 
as in the non-relativistic case: V • B = and d^/dt — V x [v x B] = (flux- freezing condition). 
Equations ((I|)-(l3|) need to be supplemented by the equation of state of the gas P = P{p,e). The 
rest of Maxwell's equations determine the current given b^ and u^: 



p-vaf) 



aOfS), 



(5) 
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w here e^'^°'^ is the Levi -Civita tensor (related to the completely antisymmetric symbol as described 



m 



Gammie et alj 120031 and references therein). Equations ([I])-® are obtained under the standard 
assumption of ideal MHD - that is, that the electric field in the fluid frame vanishes. The electric 
field in the observer's frame is non-zero, but the ideal MHD condition allows us to eliminate it from 
the equations (it satisfies E + v x B = 0), and we do not discuss it further. 

We now consider an ultrarelativistic gas in which the particle rest mass is small compared to 
the thermal energy; alternatively, our analysis is applicable to radiation-dominated flows in which 
the optical depth is high enough to ensure local thermodynamic equilibrium. The equation of state 
is then P = e/3 and the mass density can be neglected in equation In this case, equations 
([2]) and ([3]) do not contain the mass density which can be determined from the continuity equation 
given the solution u^. 

We further specialize to stationary flows {d/dt = 0). Using the condition V • B = and the 
explicit expression 7 = 1/ \/l — v^, we rewrite the linear combination of the components of equation 
([2]), -i'^v^V^T^^ - -f^v'v^V^T^^ = 0, as 

v'V'iPj'^) + -^^'V'iB"^) - j\'V^^^ = 0. (6) 
2 7^ 

In the absence of a magnetic field, equation ([6]) implies that the quantity S = Pj^ is conserved along 
any streamline and can therefore be identified with entropy, by analogy with the non-relativistic 
Bernoulli equation. This property is the consequence of our choice of the equation of state. In ad- 
dition, in the non-magnetized case the linear combination ^v'^V^T'^^ — j\7^T^^ = can be rewritten 
as V*(P^/^7f*) = 0, so according to equation ([T]) mass density in the non-magnetic fiow is oc P'^/^. 

We now consider a stationary axisymmetric [d/dip = 0) flow with a purely toroidal magnetic 
fleld [B = Bip) and a purely poloidal fluid velocity. Equations ([2]) and ([3]) can be then written in 
cylindrical coordinates z (along the axis) and r (perpendicular to the axis): 



r dr ^ ^ dz 

/.^ 9 / dvz dvz\ 9 / B^ 



2 d2\ / ^^J- dVr\ d /„ B'^\ B"^ 



0; 


(7) 


0; 


(8) 


0; 


(9) 


0. 


(10) 



^ ^ y ^ dr ~^ ^ dz J ^ dr \ ~^ 27^/ ~'~ 'y'^r 

^{VrB) + ^{VzB) 

or oz 

Equation d?]) is the 0-component of ([2]), whereas ([8]) and ([9]) are the z and r components of the 
equation y^T'^i — v^V^T^^ = and (jlOp is the (/^-component of equation 1^. The remaining 
components of equations ©-([S]) are automatically satisfied given our geometric assumptions. 

From equations ([7|)- (|10p we see that the relativistically invariant magnetization parameter 
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(the inverse of the standard ratio of gas pressure to magnetic pressure /?) quantifies the relative 
dynamical importance of the magnetic field. This value is also the ratio of the magnetic to thermal 
energy, so that Pb = 1 corresponds to energy equipartition in the rest-frame of the fluid element. 
Indeed, for a fluid element moving with 7, the rest-frame magnetic field is S/7 and the electric field 
is zero (the condition of ideal MHD), whereas pressure P is a relativistic invariant. An additional 
factor of 4tt appears in the denominator of (|lip once the magnetic field is written in Gaussian units. 
In the observer's frame, the magnetic field is B, but there is also an electric field which counteracts 
the tension of the magnetic field lines, so the dynamical importance of the magnetic tension is 
decreased by a factor of 7^. 

We now look for self-similar solutions of equations ([7|)- (|10p in the paraxial (r <C z) relativistic 
(7 3> 1) limit, and in particular we use the approximation Vz — 1 — 1/(27^) — v'^/2. For streamlines, 
the condition of self-similarity implies that fluid elements move along the curves r (x z^. Instead 
of cylindrical coordinates r and z, it becomes convenient to use the similarity variable (which is 
constant on a streamline) x = r j z^ along with z. To the lowest order in I/7, the shape of the 
streamline uniquely determines the radial velocity: 

Vr{x,z) = kxz^~^. (12) 

We then look for all solutions of the form 

7(r, z) = z^sg{x), P{r,z) = z^Pp{x), B{r,z) = z^^u{x). (13) 

After transforming equations (j7|l- (jl0p into new variables x and z (derivatives change as d/dr 
z^^d/dx, d/dz d/dz — {kx / z)d / dx) we substitute the four unknown functions in the form ()12l) - 
(fT3]) and require that x and z variables separate. To satisfy equations (f7|)- (fT0]) . the z dependence 
of the dominant terms should be the same. After this condition is satisfied, we write down the x 
component of all equations. These requirements constrain combinations of A:, /c^, kg and kp. For 
example, the presence of terms P^"^ + KB^ (where is 1/4 or 1/2) implies that j3B (defined in eq. 
[TT]) can depend only on x, so that in self-similar solutions magnetization is constant on streamlines 
(i.e., kp + 2kg = 2/c;,). We then consider three cases for Vz'- Vj- <C I/7, Vr ~ I/7 (that is, k — 1 = —kg) 
and Vr 3> 1/7. We find that there are no self-similar solutions for Vr 3> I/7, and we discuss this 
issue further in section [231 The case Vr <C I/7 yields two solutions, one with A; = and one with 
k = 1. The case t^r- ~ I/7 yields another solution with A: = 1/2. In all three cases, Pj"^ is conserved 
along the streamlines, as in an unmagnetized case. In section 12.21 we discuss the solutions with 
straight streamlines A; = (cylindrical flow) and A; = 1 (conical expansion) , and in section 12.31 the 
parabolic flow k = 1/2. 

2.2. Solutions with straight streamlines 

In the case of the cylindrical jet, k = 0, all streamlines are parallel to the axis, the radial velocity 
is Vr = 0, and the corresponding self-similar solution of (f7l)- (fT0]) is dP/dz = 0, d^/dz = and 
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dB /dz = 0. Equations ([7|)- (|10p then reduce to the weh-known (e.g., iKomissarovl Il999l ) condition 
of balance between pressure gradients and magnetic tension forces: 

d 1 d /r^^^x 

In the non-magnetic case, this condition is simply dP/dr = 0, and if it is not satisfied the radial 
pressure gradients would create radial velocities in the flow destroying the initially cylindrical flow. 
The Lorentz-factor 7(r) can vary from one streamline to another arbitrarily in the non-magnetized 
solution. 

In a magnetized flow given two of the functions B{r), 7(r), P{r) (or two independent combi- 
nations of them) we can determine the third using equation (jl4p . We assume that the magnetic 
field of the jet does not penetrate into the ambient medium (which we consider conductive enough). 
Since the total current flowing through a cross-section of radius r is /(r) = 2TrrB, this assumption 
is equivalent to requiring that the net current carried by the jet is 0. The initial condition for 
integrating equation ()14p is then 



B^ 



Pa, (15) 



where Pa is the ambient pressure. 



As an example, we consider a flow which has a constant magnetization parameter (3b{t) = const 
and is additionally homogenized so it has the same entropy on all streamlines: S{r) = P'y'^ = const. 
In this case the solution of (jl4p is formally divergent, with infinite pressure P{r) oc r"^^^/^^^^^) 
and current density j(r) = {Br)'^/r oc 7--{2+3/3B)/(2+2/3_g) ^j-^g axis, and the symmetry condition 
B{r = {)) = is not satisfied. However, the divergence is weak - the enclosed current I{r) = InrB 
and the enclosed energy flux are both zero on the axis. For a jet with a constant magnetization 
parameter, another singularity occurs at the jet boundary Vj, where a current Igc = —27rrjB(rj) 
flows with an infinite current density along the cylindrical boundary to screen the magnetic field 
B{rj) inside the jet. In Figured] we illustrate typical solutions of equation (|14p . The effect of the 
magnetic field is to pinch jet material toward the axis. 

The singularity on the axis arises because our equations (f7P- (fTU|) assume that gas pressure 
gradient is the only force that balances the pinch due to the toroidal field. Adding other pressure 
components (for example, a poloidal magnetic field) would remove the singularity. Alternatively, 
current density and pressure singularities can be removed by assuming a finite current density 
everywhere and using the current I{r) as a free function (along with another one, e.g., the same 
constant entropy condition S{r) = const as above) to determine the solution of equation (fH|) . Such 
solutions automatically satisfy the symmetry condition i?(r = 0) = 0. Since there is no current 
singularity on the boundary of the jet in this case, the initial condition ([15]) reduces to P{rj) = Pa- 
Solutions with finite current density are shown in Figure [l] with dashed lines. While they produce 
a smooth pressure profile near the axis and near the boundary of the jet, they are otherwise similar 
to the analytic solution. 
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For the solution with (3b{t) =const, the total energy flux carried by the jet is 



r°^27rrdr 



1 + ^ 27rrdr = At:^ SPa{l + PbV^, = ^T^Palh], (16) 



which needs to be multiplied by c to get the energy flux in dimensional units (cf. iLaing &: Bridle 
20021 ). Remarkably, under our assumptions (constant entropy and magnetization across the jet) 
the energy flux is independent of magnetization. For AGN jets, this gives 



140 pc X 



1/2 



5 X 10^^ erg sec~^ 



7j 



3 X 10"^'' erg cm" 



-1/2 



(17) 



The valu e for the pressure is tak en from the observations of thermal X-ray emission in the center 
of 3C31 (jHardcastle et al.ll2002l ): it is much higher than the typical interstellar pressure in the 
disk of the Milky Way (~ 5 x 10"^'^ erg cm^^) and higher yet by comparison with intergalactic 
pressures. In other words, a typ ical intergalactic niediuni cannot provide enough pressure support 
for collimation of energetic jets. iBegelman &: CiofliI (|l989l ) showed that such high pressures may be 
supplied by the cocoon which forms as the jet interacts with the ambient medium. 

Another self-similar solution of equations ([7|)- ()10p can be obtained for k = 1 and corresponds to 
a conical flow. The streamlines are straight and the variable x (which is constant on a streamline) 
is simply the angle this streamline makes with the axis, so that Vr{z,x) = x. The solution is 
P{z,x) = z~^p(x), j{z,x) = zg{x) and B = z^'^u{x). The current density is also directed along 



the streamlines, with jr{z, x) 
the conical flow is 



z ^(xu)^ and jz(z, x) = z '^{xu)'^/x. The condition for maintaining 



d 1 d 

dx^ x^ dx 



x'^v? 



0. 



(18) 



It is identical to equation ()14p . so formally the same conclusions apply for the lateral structure of 
the conical jet as for the radial structure of the cylindrical jet. However, in this geometry the jet 
is expanding fast enough that the outer parts of it are causally disconnected from the inner parts. 
We discuss the applicability of our solutions to this case in Section [27 



2.3. Parabolic solution 

There is only one more solution of equations ([7|)- (jl0p in the form (jl2p - ()13p . and it occurs at 
k = 1/2, so that the streamlines are parabolae (r = xz^/"^). In this case, the functions vary along 
the streamlines as Vr{z,x) = z^^/'^x/2^ P{z,x) = z~'^p{x), 7(2;, x) = z^^'^g{x), B{z,x) = z~^^'^u{x). 
The lateral structure of the jet satisfles 

d 1 d /x^u^\ f 2 



dx x^ dx \ 2g2 J \ 4 

Just like in the cylindrical and conical cases, the flrst two terms describe the balance between 
pressure gradients and magnetic tension forces, whereas the last term represents centrifugal forces 
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which arise because the fluid moves along curved streamlines in the parabolic case. Again, there 
are two free functions and the third can be determined from equation (|19p . In the case when the 
free functions are s(x) = pg^ = const and (3b{x) = v?' /2pg'^ = const, the solutions can be obtained 
analytically, with the initial condition (jlSp at the boundary of the jet Xj = Vj/z^/"^. Current density 
singularities similar to those discussed in the previous section are also present in this case, but since 
they are weak (in the sense that the total current and energy flowing through a finite cross-section 
are finite) and can be easily eliminated by an appropriate choice of the current density, we do not 
concern ourselves with them any further. The external pressure must decline as Pa{z) = paZ~'^ as 
the distance z from the source increases in order for the parabolic solution to be realized. Example 
solutions for the lateral structure are shown in Figure [21 

In the non-magnetized case (left panel in Figure [2]), the first free function is li = 0, and we 
again consider an isentropic flow s{x) = pg'^ = const. The solution of equation (I19p is then 



lp{x) _ ^ g]x] 



Pa 




(20) 



These solutions are shown in Figure [2] with solid lines. However, for gj > 2/xj (or alternatively 
s > Scv = ^Qpa/Xj) the right-hand side of this equation becomes negative near the axis, and 
we conclude that for these values of the Lorentz-factor (entropy) a global self-similar solution of 
MHD equations (|71)- (I10|) does not exist. The reason for this behavior is clear when we rewrite this 
condition as 'yj{z) > 2z/rj(z). In other words, the global solution in the parabolic case does not 
exist if 7 is larger than, roughly, the inverse of the opening angle of the jet 1/9 j- In this case the 
outer parts of the jet are not causally connected with the inner parts, and the lack of a global 
solution is unsurprising. 

In the magnetized case (right panel in Figure [2]) , the main difference from the purely hydro- 
dynamic case is the pinching effect of the magnetic field. The constant entropy and magnetization 
solutions are 



p{x) _ ( xy^/^'+f^^ g]x] 1 + \ g]x' 1 + /?b/2 



^Pa/{1+PB) \xjj \ 4 1 + 3Pb/2J 4 1 + 3/3b/2' 

which turn into solution (pOj) for /3b = 0. The entropy of this solution is s = Pag'j/{1 + 13b)- The 
pinching effect does not eliminate the same critical behavior as we see in the non-magnetized case, 
and for 

2 I i + Wb/2 

the global solution does not exist. 

The total energy carried by the jet can be calculated the same way as in equation (jl6p : 



Ej = ^^pagjx] - 'LpJ:±l^glx^. (23) 
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(to be multiplied by c to get the energy flux in CGS units; pa, gj and xj are all dimensional values). 
We now fix the values of Pb, Pa and Xj and consider the solutions for different values of gj. The 
behavior of the solutions is illustrated in Figure [21 with gj increasing from the top curve to the 
bottom curve in both panels. In the hydrodynamic case {Pb = 0; Figure [21 left) as we increase gj 
the total energy of the jet Ej increases all the way up to the point described by condition ([22]) . 
when the jet becomes too fast to maintain causal contact and solutions ()2ip become inapplicable. 
For the magnetized jet {(3b > 0; Figure [21 right), as we consider solutions with increasing values of 
gj, the total energy of the jet also increases at first, but then reaches its maximum value of 

Era.. = 8npa-r^^j^ (24) 

(times c in CGS units) and then decreases slightly before the jet becomes causally disconnected. 
For energies close to -Emax) most of the jet material is pushed against the wall of the jet, whereas 
for small energies most of the jet material is in the pinched core. One can think of increasing 
the Lorentz-factor of the flow as analogous to increasing the flow's inertia, so that the greater the 
inertia, the harder it is to pinch the jet toward the axis. 



2.4. Comments on flow geometry and causality 

While we do see a pinching effect of the magnetic field in the self-similar parabolic solution, 
we need to make a distinction between our solution and t he electromagnetic collimation seen in 2D 



numerical simulations. In the numerical simulations by iKomissarov et al.l (120071 ) and others the 
core is almost cylindrical when the outer streamlines are parabolic, so the ratio ^core/% decreases 
away from the base. The magnetic fields act to provide an additional collimation of the jet matter 
(the poloidal component of the magnetic field is dynamically important), not just to redistribute 
the matter close to the axis. In our case all streamlines are parabolic, and the pinched core 
also has a parabolic shape (albeit with a smaller opening angle than the outer streamlines), so 
that ^core/% =const. The difference between the geometry of our solutions and those obtained in 
numerical simulations is a direct consequence - and the limitation - of the self-similarity assumption, 
in that no additional collimation beyond what is provided by the ambient pressure can be obtained 
in the self-similar solutions. In our solution the z-dependence of all functions, including 7(x, z) = 
g{x)z^^'^ , is identical in the magnetized and purely hydrodynamic cases, indicating that acceleration 
by electromagnetic forces is matched to the thermal acceleration (their ratio is a constant on a 
streamline) . 

The lack of causal connection for > 1, which we suggested in Section 12.31 based on the 
lack of a global solution of the lateral equilibrium equation, can be understood if one considers the 
expansion of the jet cross-section in the frame moving with the Lorentz-factor 7. In this frame, 
the condition Oj^ < 1 is the same as the requirement that the cross-section expands slower than 
the speed of light. In the observer's frame, the condition 9j'j < 1 comes from requiring that the 
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time-scale for signal propagation across jet at distance z is smaller than the time scale for the fluid 
elements to get from z to 2z. 

Whether the condition < 1 is satisfied depends not just on the physical parameters such as 
energy or magnetization, but also on the general geometry. For the cylindrical solution, the value 
0j7 decreases with the distance from the origin of the jet, since 7 =const. Therefore, far enough 
from the origin the cylindrical jet is in causal contact. For the conical jet, 7 oc z and Qj =const, 
so far enough from the origin the conical jet always loses causal contact. When that happens, the 
condition Vr <^ I/7 used to derive the conical solution (end of section [2TT]) fails, since Vr = Oj. 
Therefore, the self-similar conical flow whose lateral structure satisfies equation (jlSp is no longer 
a solution of the MHD equations downstream of Ojj ~ 1. The parabolic case is unique in that 
the product Oj'j remains constant as the jet propagates. In other words, in the parabolic case the 
sound-crossing time of the jet cross-section is matched to its expansion and propagation rate. 

Our critical condition ()22p was obtained for a specific choice of free functions, namely, s{x) =const 
and I3b{x) =const. We now consider how this choice affects the presence of a critical condition 
and we take as an example a solution with a different set of free functions: g{x) =const and 
Pb{x) =const. The solution is then 

^^"^ = it:^^ UJ l2(rT^(^ - ) ■ ^''^ 

The qualitative features of this solution are the same as those of solution (|2ip. as both the edge 
pile- up (exponential term) and the magnetic pinch (power-law term) are present. Although the 
global solution formally exists for all xj and g, for large values of xjg all of the matter is piled 
up into the layer of thickness Ax ~ l/g near the edge and the pressure is exponentially small 
elsewhere. 

In this paper we investigate the structure of the jet for three special cases of the ambient 
pressure {p (x z^, p (x z~'^ and p oc z~^). In general, the interactions between the jet and 
the ambient medium may r esult in a formation of a shocked layer, such as that described by 



Bromberg &: LevinsonI (j2007l ) for a hydrodynamic stationary jet. Self-similar jets should be consid- 



ered special cases in which a solution without a shocked layer is possible. 



3. Synchrotron emission 

3.1. Constructing emission and polarization maps 

Our models relate the distribution of pressure and the magnetic field within the jet. In this 
setup, the relativistic particles that produce the pressure emit synchrotron radiation. If the total 
synchrotron losses are negligible relative to the energy flux of the jet, this emission does not affect 
the dynamical structure of the jet, and we can calculate the spatial distribution of synchrotron 
emission given P{r,z) and B{r,z). 
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We assume that all emission is optically thin and that the pressure in our model is due to 
relativistic particles with the distribution function N{E)dE = KeE~'^dE between -Emin and -Emaxj 
so the spectral index of the synchrotron emission is a = 0.5 (defined as oc and the pressure 
is P = e/3 = ln(-Emax/-E'min)/3. Since the pressure is rather insensitive to the poorly known 
value -E'maxZ-E'min! the normalization of the distribution function K^, follows the pressure during the 
adiabatic expansion as the jet propagates. Assuming a purely leptonic jet, the rest-frame emissivity 
per unit frequency and unit solid angle is 



jl,ii^',n') = 4.86 X 10"^'' |b x n'f^'^Keu'-^-^ erg s" 



^ ster ^ cm ^ Hz \ 



(26) 



where b is the re st-frame magne tic field in Gauss, is in erg cm ^ and rest-frame frequency 



is in Hz (e.g., iGinzburgj Il979l ). The unit vector in the direction of propagation of a photon 



is n in the observer's frame, but the same photon propagates with n' in the rest-frame of the 
fluid element. The rest-frame magnetic field is b = B/^ = ^/8^^J3BP, where we are now using the 
Gaussian units. For a constant magnetization parameter, the e missivity would have t raced pressure 
as ji, oc p(3+")/2 — p7/4 ^g^g ^j^g classical equipartition case, iBegelman et al.lll984l ) if not for the 
relativistic beaming and aberration effects, which are different for every volume element since they 
all move with different velocities at different angles to the line of sight. 

We integrate emissivity ()26p along the line of sight (parametrized by in the observer's frame 
and at a fixed observer's frequency v (corresponding to different frequencies z^' in the rest-frame). 
The observed flux density from a piece of jet with coordinates ^, rj as projected on the plane of the 
sky subtending a solid angle of AQ is 



(27) 



where the Doppler factor is = [7(1 — v • n/c)]~^. The aberration factor (the angle ■0' between 
b and n' in eq. 
(Appendix C of 



261) can be calculate d using the Lorentz-transformation of the propagation vector 



Lyutikov et al.ll2003l ) and turns out to be sin^ = 1 — 6"^ cos^ V- 



In addition to the intensity map ()27p . we construct maps of polarization position angle and 
polarization fraction. To this end, we calculate the emissivity into the two independent modes 
of linea r polarization, t hat is. Stokes param eters jq and jjj emitted by a unit volume within 



the jet (|Ginzburall979l : iLvutikov et al 



2OO3I ) . These expressions are similar to (j26p . but contain 



additional dep endence on the orienta tion of the magnetic field through the i?- vector of polarization. 
As detailed in iLvutikov etld] (|20oi), this latter vector also suffers relativistic aberration which we 
take into account in our calculations. The maps of (5(^, tj) and C/(^, rj) are then obtained just 
like the total intensity (127]) . Polarization fraction is \/Q(C)??)^ + U{S,,ri)'^ / S^{S,,r]) and polarization 
position angle is arctan[{7(^, r/), Q(^, r/)]/2. 

Five parameters appear in equations (j2ip . (j26p and (j27p : (j) is the angle between the jet axis 
and the line of sight; pa is the normalization of the ambient pressure; /3b is the magnetization 
parameter; Xj is the normalization of the opening angle; and gj is the normalization of the bulk 
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Lorentz- factor. Solution (j2ip has a linear dependence on pa, so changing this parameter simply 

7 /4 

results in scaling the total intensity oc pa ■ Furthermore, because of the self-similarity of our 
solutions Xj is a geometrical factor, and increasing xj results in zooming in on the region near 
the basis on the jet. Therefore, if gj is expressed as a fraction of its critical value given by (|22p 
then p{x/xj) (eq. [2T|) does not depend on Xj. As a result, only (/), 13b and gj/gj{ci) need to be 
specified to produce maps of synchrotron emission, and pa is necessary for determining the overall 
flux normalization. 



3.2. Edge-brightened jet 



As an example, we consider the jet in the nucleus of a nearby elliptical galaxy M87. The large 
mass of the supermassive black hole which powers the jet (~ 3 x 10^ Mq) and its relatively small 
distance (16 Mpc) imply that the Schwarzschild radius subtends a few micro-arcsec on the sky, 
bringing the immediate environment of the black hole within reach of modern radio observations. 
The VLBI images on scales of few h undred show that the jet is wide at the base and c ollimated 
at large distances from the nucleus (jJunor et al.lll999l : iLy et al.l 120071 : iKovalev et al.ll2007l ): in fact, 
the jet is rather close to being parabolic in shape (Figure [3^). Additionally, it is strongly edge- 
brightened, which is usually interpreted in terms of the Kelvin-Helrnholtz instability acting oii the 
boundary between the jet and the surrounding mediuin (jAloy et al.ll2000l : iPerucho et al.ll2007l : see 
also discussion of various models in lKovalev et al.ll2007l ). We propose that the edge pile-up seen in 
our models of pressure-matched jets may provide an alternative explanation for edge-brightening 
in the sub-pc jet of M87. 

Among our models, the most strongly edge-brightened one corresponds to the critical solution 
gj = gj (cr) (eq. [22l model A in Figure [2|) which we adopt as a model for the jet in M87. In this 
case the magnetization parameter f3B enters the expression for gj{cr) and therefore weakly affects 
the Doppler factor, but the pressure profile (pTI) does not explicitly depend on (3b- Therefore, 
the spatial distribution of emission is almost independent of while the total intensity scales as 
oc P'^'^ (eq. [26l) . We take Pb = 1- The last parameter needed to specify the model is the angle 
between the jet axis and the line of sight (f). In M87 it is constrained to be between 30° and 50° 
by the combination of apparent velocities of clumps (assumed to be representative of bulk fluid 



motions) and the brightness ratio of the approaching main jet to the receding counterjet (jLy et al 
20071 ). If we further assume that the jet is relativistic (v/c ~ 1) by the time it reaches the p rojected 



distance of 3.1 mas and if we take the brightness ratio measurement of lKovalev et al.l (120071 ) at face 
value {R > 200 at this distance) we obtain < 38°. We use <p = 30° in our n iodels. The observed 
emission from the jet of M87 is optically thin (a = 0.4 it 0.3; iLy et al.l 120071 ) . in agreement with 
the assumptions we made when constructing emissivity maps. 

The maps of total intensity, polarization fraction and polarization position angle for this setup 
are shown in Figure [4]A. The brighter right side is the component of the jet moving toward the 
observer. The contrast in intensity between the edge and the emission observed on-axis is similar in 
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our model and in the jet of M87 (Figure [8)3) . The brightness of the jet in the intensity map dechnes 
away from the base of the jet roughly as oc where ^ is the coordinate along the projected 

axis. This decline is due to the combined effects of the adiabatic expansion and the changes in 
aberration and the Doppler factor as the jet accelerates. The observed emiss ion declines sora ewhat 
slower, oc ^~2.4 ^ measured between the clumps from the map provided by ILv et al.l (j2007l ). One 
possible explanation for the difference in decline rates, if they are taken at face value, is that our 
models are developed in the ultrarelativistic approximation and are not strictly applicable to the 
jet in M87 which has a typical Lorentz-factor of just a few, especially near the base. Comparing 
the model intensity to the observed intensity at about 2 mas from the base of the jet (projected 
distance of 0.16 pc= SOOr^) leads to an estimate of the ambient pressure (~ 10~^ erg cm~^), of the 
total energy flux carried by the jet (a few times 10^^ erg sec~^; eq. [23]) and of the magnetic field 
(~ 1 Gauss). 

Such estimates are based on the standard assumption that the observed synchrotron emission 
from the jet reflects bulk motions and bulk pressures of the material in the jet, but it is not 
clear that this is indeed the case. An extreme example is the ra dio jet in Centaurus A, in which 
practically all of the jet emission is contained in compact clumps ([Horiuchi et al.ll2006l ). Although 
the individual exposures of the jet in M 87 are clumpier than the time-averaged composite that 
we used for our analysis (|Lv et al.l 120071 ). they all show the same edge-brightened structure in the 
process of collimation, lending some support to the bulk origin of this emission. Because the nature 
of clumps and their relationship to the bulk motions of the fluid are not understood, in our model we 
do not attempt to reproduce the som ewhat contradictory data on proper motions of these features 
(|Lv et ahlboOTl : Ixovalev et ahlboOTl ). 



Polarization measurements for the jet in M87 are now being obtained (|Lv et al.l 120071 ). and 
we present polarization maps expected from our models in Figure [HA. The polarization fraction 
of emission from every volume of fluid is the same and equals (1 + a)/(| + a) ~ 0.69, but the 
polarization position angles of emission from different elements of the volume vary, so the polar- 
ization signal washes out when magnetic fields with different directions are present along the same 
line of sight. The polarization fraction and the polarization position angle are not the same in the 
jet and in the counterjet because of the effects of aberration mentioned above, even though the 
purely geometrical projection effects are the same for the main jet and for the counterjet. In the 
edge-brightened case the polarization is expected to be high almost everywhere. The orientation of 
the magnetic field changes along every line of sight, but only a thin region near the edge contributes 
most of the emission, and the change in the direction of B across this layer is small. Stripes of 
low polarization seen in the map occur because the superposition of emission from the near and 
far walls of the jet leads to depolarization in these areas. On the projected axis of the jet, the 
magnetic field in the emitting region lies in the plane of the sky and is perpendicular to the axis, so 
the polarization E'-vector is parallel to the axis. On the edge of the jet, the projected magnetic field 
is parallel to the axis (Fig. [5]), so the polarization vector is perpendicular to the axis. In other 
words, even though the magnetic field is purely toroidal, the projection effects result in the so-called 
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'spine-sheath' polarization angle pattern (jAttridge et al.lll999l : lPushkarev et al.ll2005l ). This pattern 
is sometimes interpreted as evidence for a poloidal component of the magnetic field near the edge 
of the jet (for example, due to shearing between the jet and ambient medium). But as long as the 
jet does not lie exactly in the plane of the sky, projection effects lead to the 'spine-sheath' pattern 
even for a purely toroidal field. 



3.3. Core-dominated jet 

As another example, we now consider the other extreme in which most of the emission is 
coming from the pinched core. In this case the surface brightness on the axis formally diverges for 
> 2/5. One could eliminate this singularity by considering different profiles of Pb{x), such as 
those that we describe in Section 12. 2[ Alternatively, it can be shown that the optical depth for 
synchrotron self-absorption (opacity is oc K(,v'~^lP') diverges faster than emissivity for these values 
of (3b- So the emission from the pinched core is composed of a geometrically thin optically thick 
core surrounded by an optically thin jet. If we look at the jet at higher frequencies, we see deeper 
into the jet, so the emissivity near the optically thick core rises steeply, but at the same time the 
surface area of the optically thick core shrinks. The result is that one can always find a frequency 
at which the optically thick core contributes a negligible amount to the total emission. Observed 
jets are typically optically thin, so we assume that this condition is satisfied. The value of surface 
brightness exactly on axis is affected by the effects of optical depth, so in our calculations we assume 
that the core is smaller than a resolution element of our map and do not consider the lines of sight 
going directly through the axis. 

The resulting maps of total intensity, polarization fraction and polarization position angle are 
presented in Figure HJB. As expected, the intensity rises toward the axis. The polarization position 
angle map still displays the 'spine-sheath' structure in the main jet, albeit with a reduced 'spine', 
but in the counterjet the dominant orientation of the -E-vector is now perpendicular to the axis. 
The 'spine' shrinks because the emissivity rises toward the axis, so the largest contribution to the 
surface brightness is acquired where the line of sight passes closest to the axis, where the magnetic 
field projected on the plane of the sky is entirely parallel to the axis (tangential points in Figure 
[5]). The naive conversion of the polarization position angle into the orientation of the magnetic 
field would imply that the field is mostly parallel to the axis, but in fact the fields are purely 
toroidal. This example illustrates that not only the orientation of the jet is important in inferring 
the direction of the magnetic field, but the emissivity profile as well. 

In the core-dominated case, the polarization fraction is small in the center, where multiple field 
directions contribute, resulting in depolarization. The polarization fraction is close to its theoretical 
maximum on the edge, where only a thin layer of the jet material contributes to the observed 
emission. Thi s rise is reminiscent of the puzzling pola rization fraction profil es in the parsec-scale 



jets of 3C273 (IZavala fc Tavlod I2nn5l ) and in Mkn501 (|Pushkarev et al.ll2nnRl ). We suggest that in 



these objects the emissivity is rising toward the axis, producing the observed polarization pattern. 
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In both these objects the orientation of the i?-vector of polarization is perpendicular to the jet in 
the regions where the emission is highly polarized, in agreement with our predictions. To further 
strengthen the case for the toroidal fields in these objects, Mkn501 shows a well-defined 'spine- 
sheath' polarization pattern, while 3C273 shows strong gradients of the rotation measure across 
the jet. 

The core-dominated case is only mildly relativistic (the highest Lorentz-factor on the edge of the 
jet is only 7 ~ 3 in FiguredP), so our solutions for the lateral equilibrium are not, strictly speaking, 
very accurate in this case. However, the effects of the geometric projection and emissivity profile 
on emission and polarization maps are rather generic and are well-illustrated by this approximate 
model. 



4. Conclusions 

We have considered axisymmetric relativistic MHD jets with an ultrarelativistic equation of 
state and toroidal magnetic fields. We have then found all self-similar solutions of relativistic 
MHD equations in this setup assuming that the jet propagates in a pressure equilibrium with the 
surrounding matter. There are three geometrically self-similar solutions which occur in the medium 
with pressure scaling as P{z) oc (cylindrical jet), oc z~'^ (parabolic jet) and oc z"^ (conical jet). 
The lateral structure of the jet is determined by the balance of the radial pressure gradient, the 
tension of the toroidal magnetic field lines and the centrifugal force. The last term is relevant only 
when the streamlines are curved. Our models represent a special case of a m ore general situation in 



which interactions between the jet and the medium lead to shocked layers (iBromberg &: Levinson 



20071 ): in our case, solutions with no shocks are possible. 



While the self-similarity condition is quite restrictive, in that only a limited set of geometries 
can be studied, the lateral structure of the parabolic jet displays some generic features which we 
investigate in detail. The pressure gradients act to push the jet material against the walls of the 
jet, whereas the magnetic field tension acts to pile the material closer to the axis. The relative 
importance of the two effects is determined by the physical parameters of the jet, such as its 
magnetization Pb and the total energy. The pinched core component is stronger for a higher value 
of Pb, whereas the pile-up component is stronger for a jet with a higher Lorentz-factor (or higher 
total energy, to put it another way). 

Our models describe jets whose Lorentz-factors are smaller than the critical value 7cr — 
where 6j is the opening angle of the jet. We find that for larger Lorentz-factors 7 ^ 7cr either 
the formal solution of the lateral equilibrium equation does not exist at all, or it exists in the form 
of a thin pile-up of the material along the edge of the jet. The condition 7 < 7cr is equivalent to 
requiring that the jet maintains causal contact in the lateral direction as it propagates. 



Edge pile- up is a generic feature of a jet with curve d streamlines and it might explain edge- 
brightening of the jet in M87 in the acceleration region (jLv et al.l 120071 : iKovalev et al.l 120071 ). We 
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also find that even purely toro idal magnetic fields can produce the 'spine-sheath' polarization 
structure (jPushkarev et al.ll2005l ) when projection effects are taken into account. We find that for 
core-dominated jets the polarization fraction is expected to be near ze ro on the projected jet axis 
and rise toward the ed ge, as observed in 3C273 (jZavala Tavloilbood l and Mkn501 (B1652+398; 



Pushkarev et al.ll2005l ). The orientation of the polarization position angle is determined not just 



by the geometric projection effects, but also by the emissivity profile. For emissivity rising toward 
the axis, the polarization perpendicular to the axis may become the dominant orientation even 
for purely toroidal magnetic fields. These effects significantly complicate the interpretation of the 
observed polarization position angle in terms of the intrinsic direction of the magnetic field. 
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Fig. 1. — Lateral structure of a cylindrical jet: non-magnetized (dotted line), constant magnetiza- 
tion parameter Pb (solid line), /3b (r) chosen so that the current density is finite everywhere (dashed 
line). Left panels show a weakly magnetized case {Pb ^ 0.1) and right panels show a magnetized 
case {Pb ^ 1.0). The weakly magnetized solutions differ from the magnetized solutions only by a 
constant stretch since for a given /3b (r) the equation to determine pressure (jlSp is linear. The jet 
boundary is at rj = 5. 
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Fig. 2. — Lateral structure of a parabolic jet: non-magnetized (left) and magnetized with a constant 
magnetization parameter [3b = 1 (right). In each panel, the different lines, from top to bottom, are 
analytical solutions of equation (|19p with increasing value of the Lorentz-factor (parametrized by 
Qj), with all other parameters {(3b ^ Pa and Xj) fixed. The dotted lines show the cylindrical solution 
for comparison; the solid and dot-dashed lines show physically acceptable solutions gj < gj{cr), 
where the critical value is given by (j22p . The dashed lines show an unphysical solution gj > gj{cx) 
for which pressure formally becomes near the axis. The thick solid lines show the critical solution 
with gj = 5j(cr). In the non-magnetized case (left) the energy carried by the jet monotonically 
increases downward, whereas in the magnetized case (right) it reaches its maximum value given by 
equation (j24p on the dot-dashed curve and then declines. Labels A and B indicate the models used 
in calculation of synchrotron emission in Figure HI 
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Fig. 3. — (a) Opening an gle of the jet in M8 7 as a function of distance from the core. The points 
with error bars are from lJunor et al.l (119991 ) and references therein. For comparison, the dotted 
hne shows the opening angle of a parabolic jet {0j oc l/-v/i); it is not a fit to the data, (b) Latera l 
surface brightness profile of the jet in M87. The profile constructed from data by lLy et al.l (120071 ) 
is shown as a histogram. The profile from our model is shown with a thick solid line. Effects of the 
finite beam size have not been taken into account. 
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Fig. 4. — Maps of synchrotron emission from our models. Model A is for the critical solution 
(thick solid line in Fig. [2)3), in which all of the material is piled up against the wall of the jet, 
whereas model B is for the core-dominated solution (top solid line in Fig. [2)3). For each of the two 
models, we show intensity (square root scaling), polarization fraction (blue - low, red - high) and 
polarization position angle counted clockwise from the vertical direction (yellow - polarization is 
along the axis, blue - polarization is perpendicular to the axis). The core region is excluded from 
the maps. The right side of the jet is moving toward the observer. The Lorentz-factor of the fluid 
follows 7 oc P-I/^; the normalization is such that for model A, 7 ~ 15 for the fluid elements in the 
top right corner of the intensity map, and for model B, 7 ~ 3 for the same elements. 
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Fig. 5. — Appearance of a toroidal field line (thick) and the axis of the jet (thin) as projected on 
the sky for a set of inclination angles ^. Even though the field is purely toroidal, for all angles 
(f) 7^ 90° it will appear as parallel to the axis in the tangential points Ti and T2. 



